Generic evaluation of the relaxation time to equilibrium 
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We evaluate the relaxation time to equilibrium, and especially show that it is almost independent 
from the system size for macroscopic isolated quantum systems. It at most polynomially depends on 
the system size. This estimation holds when the Hamiltonian is non-integrable, the initial deviation 
of the quantity of interest is of order its spectral norm, and the relaxation process is monotonic. 
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Recently many attentions have been paid for the foun- 
dation of equilibrium statistical mechanics. For macro- 
scopic quantum systems, the relaxation to equilibrium 
has been explained by probabilistic arguments based on 
the extreme high dimensionality of the Hilbert space 
stems from the noncommutativitvfll-fl2l) . Refs-fH-Q es- 
tablish that majority of the states are regarded as equi- 
librium, whose properties are commonly shared. On the 
other hand, Refs.[y? 413 ] explain the relaxation process 
by also assuming either macroscopic properties of observ- 
ables or nonintegrability of the Hamiltonian. 

But the evaluation of the relaxation time remains as a 
significant unsolved problem. It would strongly depend 
on models and quantities of interest, and thus is believed 
to be difficult to estimate. Even very fundamental prop- 
erties are in general unclear such as how short the relax- 
ation time is compared with an extremely long recurrence 
time[13j, or how strongly it depends on the quantities of 
interest given that the relaxation occurs. 

In this letter, we show a generic estimation of the re- 
laxation time under the main assumptions i-iii). The 
assumptions are as follows. 

i) The Hamiltonian is nonintegrable. 

ii) For a quantity of interest A which polynomially 
depends on the system size, the deviation 5A(0) 
of the initial expectation value from the long time 
average is less than a kind of spectral norm ||^4||. 

iii) In the course of time evolution, the deviation of the 
expectation value SA(t) monotonically approaches 
to zero. 

The assumption iii) implicitly assumes that the relax- 
ation occurs. Physically, the relaxation time is expected 
to be almost independent from the system size or vol- 
ume for macroscopic systems. Indeed, macroscopic sys- 
tems are divided into many pieces which are mutually 
interacting, and each portion would relax independently. 
From this macroscopic argument, the polynomial depen- 
dence at most would be trivial, however, here it is shown 
as a generic theorem. 



Let us consider a normalized state \tp{t)) in a large 
Hilbert space H[e,e+ae] a t an energy E. We are con- 
cern with the relaxation time for a large but isolated 
system. Justifications for the use of isolated system to 
describe relaxation has been discussed for example in 
Ref . fiol] . From an initial nonequilibrium, it evolves by 
a time-independent Hamiltonian H whose eigenenergies 
and eigenstates are E n and \E n ), respectively 
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where c„ and <p n are the amplitude and phase of eigen- 
state \E n ). d = <&bM\e,e+&.E\ 1S the dimension of the 
Hilbert space, which exponentially depends on the sys- 
tem size N. The width AE is smaller than E, but much 
larger than neighboring energy distances of order ^jp. 
Such a width is always necessary for the microcanoni- 
cal approach. We consider a physical quantity A whose 
maximum value \\A\\ = Max|$) eW[E B+AE] polyno- 
mially depends on the system size N. For example, the 
sum of local finite quantities is proportional to N and sat- 
isfy this condition, and these are important macroscopic 
quantities \§L fl2j. 

The difference between the expectation value and long 
time average of a quantity of interest A is calculated as 

SA(t) 
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where we evaluated c n = 0(^=) from the normal- 
ization, and |(I£ n |j4|i? m )| = 0(^=1) for nonintegrable 
systems [12[, which is correct including the standard devi- 
ation. Regarding the evaluation of c„, it is also possible 
to say that we choose c„ = and rigorously evaluate 
the order of the relaxation time. The broad distribution 
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over the eigenstates is also reasonable for the realistic 
preparation of the initial states[l(|. 7 n , TO = — 7 m ,n is the 
phase of the off-diagonal element (E n \A\E m ), which is 
in general a function of (m, n). Note that the diagonal 
contribution expresses the long time average 
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Also note that our analysis only requires that there is a 
typical small amplitude of off-diagonal elements. And we 
shall skip the quantitative evaluation of the off-diagonal 
elements since it is detailed in Rcf. 12] with numerical 
evidences, but its smallness is explained as below. In- 
tuitively, the off-diagonal elements are extremely small 
compared to the diagonal ones, since they are related to 
the transition amplitudes between macroscopically dif- 
ferent states \E n ) and \E m ) by a perturbation A, which 
is an almost impossible process. Note that the argument 
here holds also for E n very near to E m , since the shape 
of the corresponding eigenvectors are different. 

The initial state shows deviation from equilibrium 
value, and we assume 



For a — | — e, there is a strong correlation among the 
phases {4> m } compared with uncorrelated case a = 1. 
In this paragraph, we estimate how long the phase cor- 
relation is maintained. Phases S — — 4>n ~ 'Pm + 
1mm < 777, 77 < d} are decomposed to Si whose elements 
are concentrated around a specific value < 4> < 2tt and 
a completely uncorrelated set S2 = S—Si whose elements 
randomly distributed over [0, 27r]. Strictly speaking, 
there can be several <fr, around which the correlated com- 
ponents of Si distribute. But for simplicity, we assume 
single valued (f>, which does not affect our analysis. In the 
course of time evolution, Si is well-defined when the ad- 

(3) 

ditional phases satisfy I- 5 " n Ern ^ t <C (j) for most 77 and 777. 
This condition is satisfied for t < -Pw(f>, t ne decorrelation 
time scale for 77 and 777 with \E n — i? m |near AE. Other 
additional phases with much smaller \E n — E m \ < AE 
does not affect the correlation of Si within the duration. 
We also required that £2 remains uncorrelated. 

Let us calculate the relaxation time T r . The time 
derivative of the expectation value is calculated as 
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where the important parameter a is determined later. 
Note that if the phases {<p n — 4> m + 7, lm } are uncorre- 
lated, the sum over the random phases corresponds to 
d(d— 1) steps of a unbiased random walk, and thus a = 1 . 
To clarify this point, let us consider a uniformly random 
quantity £„ m and estimate | , e^™ m |. For a fixed 

777 = 7770, Ree l ^" m o is also unbiased random variable tak- 
ing positive and negative values with an equal probabil- 
ity. The summation over 77 is therefore of order yd and 
can be both positive and negative. Then the sum over 
777 is actually regarded as a random walk with an almost 
common step Vd, which yields total displacement O(d), 
i.e. a = 1. 

Substituting the ansatz Eq. (HJ) to the evaluation of ex- 
pectation value Eq.([2]), we have 



|(^(0)|A|^(0)} - £ cl(E n \A\E n )\ = d«-i\\A\\. (5) 

71 = 1 

We assume that the initial deviation 5A(0) is of the same 
order as \\A\\. Therefore the exponent satisfies a — | — e 
with a very small positive quantity e = +0. e accounts 
for the case that the deviation is slightly smaller than 
\\A\\. We use this notation for e, since d — 0(e N ) is 
extremely large number and e can be smaller than any 
positive number for sufficiently large N. 



f° r t _• ~Ke4>- This provides the rate of relaxation. Now 
we assume that the expectation value monotonically ap- 
proaches to the equilibrium. The relaxation time re- 
quired for this change of the expectation value is eval- 



uated as T, 
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= j \A j j from Eq.©, and thus 



Tr = ^d*-*. (7) 

The relaxation time is therefore evaluated as 



Tr = (8) 

Therefore from the definition of an extremely small pa- 
rameter e, the dependence on d of the relaxation time is 
weaker than any polynomials of d. The logarithmic de- 
pendence on d is possible, and T r at most poly normally 
depends on the system size N. This includes the real- 
istic case where the relaxation time T r is independent 
from system size. More generally, T r = ^N a with a > 



if we choose e = 



a log 1 



which can be smaller than any 



logd 

positive quantity as N — > 00. 

In conclusion, we give the first step of the studies for 
generic properties of the relaxation time. We evaluated 
the relaxation time and showed that it at most polyno- 
mially depends on the system size for generic systems 
satisfying assumptions i-iii). The evaluation guarantees 
a reasonable relaxation time without going into the detail 
of a specific model. 
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